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Theoretical investigation of the dynamic electronic response of a quantum dot driven 

by time— dependent voltage 
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(Dated: submitted on 22 July 2008s) 

We present a comprehensive theoretical investigation on the dynamic electronic response of a 
noninteracting quantum dot system to various forms of time-dependent voltage applied to the 
single contact lead. Numerical simulations are carried out by implementing a recently developed 
hierarchical equations of motion formalism [J. Chem. Phys. 128, 234703 (2008)], which is formally 
exact for a fermionic system interacting with grand canonical fermionic reservoirs, in the presence of 
arbitrary time-dependent applied chemical potentials. The dynamical characteristics of the transient 
transport current evaluated in both linear and nonlinear response regimes are analyzed, and the 
equivalent classic circuit corresponding to the coupled dot-lead system is also discussed. 

PACS numbers: 05.30.-d, 72.10.Bg, 73.63.-b 



I. INTRODUCTION 

Motivated by the rapid development in the field of na- 
noelectronics, a comprehensive and fundamental under- 
standing of quantum transport phenomena has become 
an urging quest. Theoretical investigations on electronic 
dynamics of open quantum systems subject to external 
fields not only provide great insights into relevant physi- 
cal problems, but also shed light on the design and ma- 
nipulation of mesoscopic or nanoscopic electronic devices. 

Frequency-dependent admittance of mesoscopic sys- 
tems has been studied by scattering theorji^i^ as well 
as nonequilibrium Green's function (NEGF) methodi^i^i^ 
At low frequency and low temperature, it was found that 
to linear order the coherent quantum dynamics of a meso- 
scopic capacitor can be characterized by a classical cir- 
cuit in which a resistor and a capacitor are connected in 
series. '^'^'■^''^ The resistor Rq, often termed as the "charge 
relaxation resistor" , is related to the time scale for elec- 
trons to accumulate on the capacitor in response to ex- 
ternal applied voltages, and does not rely on transmission 
coefficients. For a single-channel, spin-polarized contact 
with a single lead, it was predicted that at zero temper- 
ature Rq is equal to half a resistance quantumji i.e.. 
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(1) 



The factor 1/2 arises due to the fact that there is only one 
channel connecting the quantum system to single elec- 
trode. This has been confirmed quantitatively by a re- 
cent experiment.— The capacitance C^, normally referred 
to as the "electrochemical capacitance" , was found signif- 
icantly different from the electrostatic geometric capaci- 
tance Co due to the limited size of the quantum system 
and the coherence nature of electron transport. The clas- 
sical circuit was extended afterwards to include an addi- 
tional "quantum inductance" Lq, if one of the capacitor 
plates is a quantum dot (QD)i^i^i^ i-^°i^^ Lq is of purely 
quantum origin, and is associated with the time scale of 
the resonance. For instance, consider a QD consisting of 
a single spin state and coupled to one electrode. Its dy- 
namic admittance at zero temperature, G{uj), has been 
expressed analytically by mapping to the equivalent cir- 
cuit with resistance Rq, capacitance C^, and inductance 



l^q a^" 

G{uj) = -iujC^ + uj^C^Rq + iLo'^ClRl - iuj^ClLq. (2) 

So far, most work in this field has been conducted 
within the linear response regime and based on analysis 
in frequency domain. Keeping in mind that the electron 
transport is actually a coherent process taking place in 
real time, it is thus intuitive and straightforward to study 
the dynamic properties of an open quantum system by 
probing its transient response to external applied fields. 
The frequency-dependent admittance G{uj) can be ob- 
tained via 



G(c.) = 



/(c) ^ T[m] 

Vic) :F[v{t)r 



(3) 



where J^ denotes conventional Fourier transform. It is 
worth emphasizing here that G(w) is actually indepen- 
dent of the applied voltage V{t), so long as its amplitude 
is kept sufficiently small. Equation ([3]) shows an alter- 
native route to evaluate G(w) provided that the time- 
dependent current response I{t) can be accurately simu- 
lated via for example quantum dissipation theory (QDT). 

We have recently constructed a formally exact QDT, 
in terms of hierarchical equations of motion (HEOM), for 
arbitrary non-Markovian dissipation systems interact- 
ing with Gaussian grand canonical ensembles i^^'^'^i^^'^^ 
The theoretical construction was carried out on the ba- 
sis of the calculus-on-path-integral algorithm, together 
with the spectral decomposition technique] ^'^i^^'^^i^^'^^ 
Dynamic responses of the reduced quantum system to 
external fields can then be obtained through numeri- 
cal solutions of HEOM for the system density matrix 
and its associated auxiliary counterparts. For a general 
many-particle system, the HEOM formalism needs to 
be properly truncated at a certain finite tier, and the 
simulated outcomes are considered reliable as long as 
they are convergent with respect to further inclusion of 
higher tiers F^^ii^ii^ii^ Great simplification does exist for 
single-particle systems. It has been proved that for the 
electron transport through a noninteracting system, the 
HEOM formalism terminates strictly at the second tier 
("max = 2) without approximation 3 

In this work, we focus on the electron dynamics of a 
single-level noninteracting QD coupled to an electrode 



with finite bandwidtli. The paper is organized as fol- 
lows. In SecHH the theoretical framework and prac- 
tical implementation of the HEOM formalism is intro- 
duced. In Appendix, numerical accuracy of the HEOM 
approach is validated by comparing to exact quantum 
transport calculations reported in literature. In Sec. lIIII 
the HEOM approach is applied to simulate the transient 
electronic dynamics of an open QD. Numerical results 
in linear-response regime are presented and discussed, 
along with a detailed analysis of frequency-dependent 
dynamic admittance. In Sec. IIVI nonlinear effects are ex- 
plored. Transient current responses to various types of 
applied voltages will be exemplified. Conclusions and 
further comments are given in Sec.FVl 



II. METHODOLOGY 

The QDT-HEOM formalism developed recently is for- 
mally exact for the following standard form of Hamilto- 
nian for quantum transport,— 



H^ 



H{t;{a^s,aU)+y^{h^+H'J. 



(4) 



The electronic Hamiltonian (H) of the system (such as 
QDs) is rather general, including Coulomb interaction 
and time-dependent external fields. The electrodes are 
modeled by noninteracting electrons, 



'T'a — / ^ / ^ ^aks^Q^j^g^aks- 



(5) 
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The transfer coupling i/^ between the system and the 
a-electrode reads 



kea fJ.s 



<ks"- 



flS 



H.c. 



(6) 



Here, afj,s (atg) is the annihilation (creation) operator 
associated with the single-electron state /x and spin s 
of the system, daks {d-aks) ^^ ^'^^^ associated with the 
specified a-electrode single-electron state of energy Caks , 
and takfis is the transfer coupling matrix element. The 
transfer coupling spectral density functions are 
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(7) 



In contact of QDT, which describes the dynamics of re- 
duced system density operator p{t) = tisPTit), we treat 
the electrodes as the grand canonical fcrmionic reser- 
voir bath. We denote hs — J^a ^a and the thermo- 
dynamic density operator p'^ for the bare bath in the 
absence of time-dependent bias voltage. To describe 
the stochastic nature of the system-reservoir coupling, 
consider Eq. ([S]) in the /iB-interaction picture, H'^{t) — 
E^s flfisitWs + H.c, with (setting ft = 1) 



/t^,(t)^e'"B*(^^^^ 
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These stochastic bath operators satisfy the Gaussian 
statistics with Wick's theorem for thermodynamic av- 
erage over the grand canonical fermionic bath ensem- 
bles with p'^. The effects of electrodes on the system 



are completely the following nonzero reservoir correlation 
functionSfi^ 

CU.sd -r)^ {fl^sit)fa.s{T)U, (9a) 

C-^.sit -r)^ {Us{t)fLsir)U, (9b) 



which relate to the spectral density functions (denoting 

Jafj.iys) Via — 



7+ = J- 



i noo ±iuj{t-T) j± ( \ 
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at,,ys\ ) 2tT J_^ l+e±/3„(^-M=) ^ ' 

Here /?„ = l/{kjiToi) and pa are inverse temperature and 
equilibrium chemical potential of lead a, respectively. 

In the presence of time-dependent voltage applied to 
electrodes, the Fermi energy is subject to a homogenous 
time-dependent shift, i.e., Pa{t) = Pa + Aa(t). The 
relevant nonequilibrium correlation functions arei^ 



Canvsit^'^) =exp 



±i / dt'Aait') 



Ct,.sit-r). (11) 



From the perspective of QDT, the electronic dynamics 
of an open system is mainly characterized by p{t) and the 
transient current through each lead. For linear coupling 
Hamiltonian of Eq. ([6]), system-lead dissipative interac- 
tions can be exactly captured by the Feynman- Vernon 
influence functional in path integral. ^^ Time derivatives 
on the influence functional are performed in a hierarchical 
manner, and thus lead to the construction of a formally 
exact hierarchical set of coupled EOM for a general non- 
Markovian dissipative system F^^ii^iiiii^ The final HEOM 
is cast into the following compact form (see Ref.[l3J for 
details) : 



Pn = -[iC + 7n(<)] Pn+ Pn + Pn 



,{ + } 



(12) 



with CO = [H{t),0] for an arbitrary operator O. The 
basic variables of Eq. (fT2|) are p{t) and associated auxil- 
iary density operators (ADOs) Pn{t), where n is an in- 
dex set covering all accessible derivatives of the influence 
functional. With C^ausi^) expanded by an exponential 
series via a spectral decomposition technique^^iiiii^ n 
involves an n-fold combination of {cr,a, p,i>, s,m) that 
characterize the exponential series expansion with cr — ±. 
Therefore, pnU=o = p{t), and 7n|fi=o = Pn"Ms=o = 0; 
Pn\n^Q is an ADO at the n^^-tier; 7n(i) collects the re- 
lated exponents along with iAa{t), to p^; and pn~^ and 
Pn^' are the nearest lower- and upper-tier counterparts 
of Pn, respectively. In particular, the I'^^-tier ADOs, 
Pn(OU=i = Pa/ii^smi^) ' determine exclusively the tran- 
sient current of spin-s through the lead a asF^ 



/„,(f) = -2Im^tr[ 



1+ 



{t)] 



(13) 



Here the trace is performed for all QD degrees of freedom. 

The HEOM formalism has been implemented for a gen- 
eral dissipative system. Details of the programming tech- 
niques will be published elsewhere. The accuracy of our 
numerical approach is verified extensively; see Appendix. 

In the following we will focus on a single-level non- 
interacting QD coupled to one lead, and thus omit the 



indexes of (a^vs) hereafter. The energy of the spinless 
QD-level is eo- An external voltage V{t) is applied to the 
coupling lead from i = 0, which excites the QD out of 
equilibrium. The lead energy level are shifted due to the 
voltage, A(i) — —eV{t), and so is the lead chemical po- 
tential fi{t) = fi + A(t), where e is the elementary charge 
and fi is the equilibrium lead Fermi energy that is set to 
zero hereafter, i.e., /i = 0. A widely used Drude model is 
adopted for the spectral density function of the coupling 
lead, 



J{u;) 
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In equilibrium (in absence of voltages) or in a steady state 
(under a constant external voltage), the rhs of Eq. (|12|) 
is equal to zero, and the HEOM reduce to a closed set of 
linearly coupled equations for {pn{0)}. The subsequent 
evolution of {p„{t > 0)} is characterized by Eq. P^ . 
with the equilibrium density operators as initial condi- 
tions, i.e., Pn(0) = Pn'^. With our current coding, the 
reduced system is spanned in a Fock-state representa- 
tion. The linear sparse problem for solving the equilib- 
rium or stead-state pn is tackled by the biconjugate gra- 
dient method)^ and the time propagation of p„ (t) in fol- 
lows the 4"^-order Runge-Kutta algorithm. It has been 
shownH, that for an noninteracting system, the HEOM 
formalism is exact with the terminal tier of nmax — 2. 



III. DYNAMIC ADMITTANCE IN LINEAR 
RESPONSE REGIME 
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FIG. 1: Transient current responses to various types of exter- 
nal voltages. The solid (dashed) curves represent the time- 
dependent current I{t) [energy shift A(t) = —eV{t)]. The pa- 
rameters adopted are as follows (in unit of meV): eo = /^ = 0, 
r = 0.08, r = 0.026 and W = 3. The switch-on vohages 
are (a) a step function; (b) a square pulse; (c) an exponential 
function; (d) a sinusoidal function; (e) a symmetric Gaussian 
function; and (f) an asymmetric Gaussian function with the 
same peak amplitude A — 0.001 meV. 



A. Transient electronic dynamics and 
frequency— dependent admittance 

We now apply the HEOM approach to investigate the 
dynamic admittance of the single-lead QD. It is pre- 
sumed that the geometric capacitance can be neglected 
for the QD of interest, i.e., Cq — > cx), which means that 
no charging energy is required for an electron to popu- 
late the QD level. It is also assumed that the QD level 
energy eo{t) = eo -|- Aoit) does not change with time, 
i.e., A]j{t) = at any time t, which is consistent with 
the Co — > oo hypothesis. The dynamic admittance of a 
single-lead QD is extracted from the real-time current 
responses I{t) to external driving voltages V{t) based on 
Eq.®. 

The general HEOM formalism admits an arbitrary 
form for the external voltage V{t). Transient currents 
due to various types of applied ac voltages are shown in 
Fig-El The time-dependent voltage adopted for Fig.[IJc) 
corresponds to A{t) = A(l — e~*/'^"), where the time 
constant Tq > dominates the switch-on rate of the 
voltage, and the asymptotic limit Ta -^ 0"*" actually de- 
scribes a step function. In Fig.[IJc) Tq is taken as 50 ps. 
The driving voltage for Fig.[TJd) is a sinusoidal function: 
A{t) = A[l - cos(a;i)] with the period (27r/cj) = 40 ps. 
In Fig.[lje) and (f) the open QD system is excited out of 
equilibrium by a Gaussian voltage pulse as follows. 



A(t) 






<t <Tc 
t> Tc 



(15) 



where r;, and Tc determine the width and the center of the 
Gaussian pulse, respectively, and the factor k controls the 
asymmetry of the pulse before and after Tc. In Fig.[lje) 
Tb = 22 ps, Tc = 100 ps and k = 1 are employed, i.e., the 
Gaussian pulse is symmetric in time, while in Fig.[TJf) 
K = 4 is adopted, which means that A(t) drops faster 
than it rises. In the linear response regime, the dynamic 
admittance of the open QD system obtained by Eq. ^ 
reflects the intrinsic physical features of the open system, 
and is thus independent of the specific temporal behav- 
ior of A(i), as long as its amplitude A is kept sufficiently 
low. Among the various time-dependent voltages men- 
tioned above, the asymmetric Gaussian pulse is found 
to be a convenient candidate for a Fourier analysis, since 
both V{t) and the corresponding I{t) have nonzero values 
only within a finite time interval. Therefore, subsequent 
numerical analyses in this section will be based on sim- 
ulations with asymmetric Gaussian voltage pulses, if no 
additional remarks are given. The purpose of introducing 
asymmetry to the Gaussian voltage pulse [see Eq. (fT5|) ] is 
to save computational effort while maintaining the accu- 
racy of the resulting admittance G{uj), especially in the 
low frequency range (to be elaborated later). Transient 
current driven by some other types of applied voltages 
will be discussed in Sec. lIVi such as step-function and 
delta-function voltage pulses. 

It is noted that for a voltage pulse consisting of both 
turn-on and turn-off sides, a sign change is always ob- 
served for the response current [of. Fig.[T]^b), (e) and (f)], 
which shows that an entirely positive A(t) can result in 
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FIG. 2: Transient current responses to an asymmetric Gaus- 
sian voltage pulse (see the inset) under various temperatures. 
The parameters are « = 4, and the rests in unit of meV: 
A = 10"*, r = 0.1 and W = 10. 
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a negative current. This is actually due to the fact that 
there is only one lead coupled to the QD. As the voltage is 
turned on, the lead chemical potential is increased, which 
drives the electrons flowing from the lead to the QD (pos- 
itive current) until the voltage reaches its maximum; and 
then while the voltage decays, the excess electrons resid- 
ing on the QD gradually wane back to the lead, which 
reverses the direction of the electron flow (negative cur- 
rent). As the voltage pulse vanishes at i — > cx3, the initial 
equilibrium is restored ultimately. Therefore, during the 
entire voltage on-off cycle there is no net charge accumu- 
lating on the QD. Based on the charge continuity equa- 
tion, we should have I{lo = 0) = /_ /(t) dt = 0, and 
this equality has been verified for all the cases plotted 
in Fig.[Tl except for (d) where A(t — > oo) ^ 0. For the 
case of an asymmetric Gaussian voltage, it is intriguing 
to see the negative peak current (31.5 pA) is larger than 
the positive counterpart (25.8 pA); see Fig.[Hf). This is 
ascribed to the non-adiabatic charging effect. In gen- 
eral, the more rapidly the external voltage changes, the 
larger is the transient response current in terms of its 
amplitude. Here the turn-off side of voltage possesses a 
steeper slope (k > 1), and hence the electrons going out 
from the QD at i > Tc is faster than the rate of electron 
inflow at i < Tc. For the same reason, the peak current 
under a step function voltage is larger than that under 
an exponential function voltage; c/. Fig.[lja) and (c). 

In linear-response regime, the dynamic admittance of 
a noninteracting QD coupled to a single lead has been 
derived by the NEGF method with the wide-band limit 
(WBL) approximation.^'^ Extension to a finite band- 
width case is straightforward. With a Lorcntzian spectral 
density function Jit) [c/. Eq. ((M])], the linear-response 
admittance G{lS) at any temperature can be evaluated 
via Eq. ^ for both resonant (eo — /i) and off-resonant 
(eo 7^ /i) cases with /i being the equilibrium lead chemical 
potential. 



FIG. 3; (a) Real and (b) imaginary parts of frequency- 
dependent admittance under various temperatures. Same pa- 
rameters are adopted as in Fig. [21 



B. Resonant tunneling cases 

We first investigate the cases where the dot level eo 
is in resonance with the Fermi energy of the lead /i, i.e., 
eo = /J. = 0. In Fig.[2]we plot the transient currents driven 
by an asymmetric Gaussian voltage pulse under various 
temperatures. For all calculations carried out, the ampli- 
tude of applied voltage is kept lower than 10 /iV through- 
out the simulation time to ensure the system dynamics 
remains in the linear response regime. Upon the applica- 
tion of the driving voltage, the QD undergoes a complete 
period of charge accumulation and depletion, which is 
synchronized with the on-off cycle of the external bias. 
It is observed that the current amplitude is suppressed 
as the temperature rises. The corresponding frequency- 
dependent admittances are shown in Fig. [31 where the 
conductance quantum 50 — 'i-e^/h = 7.75 x 10^^ S is 
used as the unit for G{ljj). It is noted that the imaginary 
part of G(w) changes its sign within the frequency range 
uj/{2i:) G (30, 40) GHz [see Fig.lSIa)]. This implies that 
the electronic dynamics of the reduced system exhibits 
distinct phase features at different frequencies. In the 
low frequency regime, Im[G(w)] is almost proportional 
to w, hence the open QD system resembles much a ca- 
pacitor. Since the geometric capacitance Co is neglected 
for the open QD system, this capacitor-like behavior is 
of pure quantum coherence nature, and is thus referred 
to as electrochemical capacitance C^. As a; increases, the 
linear relation for Im[G'(a;)] breaks down, and the phase 
shift between the current and voltage diminishes grad- 
ually until getting reversed upon the sign change. This 
agrees with previous studies on the C^ of the same open 
QD system»^ 
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FIG. 4: 'Re[Z{uj —^ 0^)] as a function of F under various 
temperatures. Other parameters are (in unit of meV): A = 
10~^ and W = 10. The inset plots Rq as a function of a 
dimensionless quantity F/T. 
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FIG. 5: Re[Z((i;)] as a function of the dimensionless quantity 
ui/r. The lines represent the different values of F. Other 
parameters are (in unit of meV); T = 0.026 and W — 10. 



At zero temperature, it has been demonstrated that 
the dynamic admittance of the open QD system can be 
quantified by a classical RLC circuit in the low frequency 
range (w < F)!^ Adopting the convention in quantum 
transport theory, the electrical impedance of the RLC 
circuit is expressed as follows [cf. Eq. ^], 



Z{uj) = [G{uj)]-^=R,+ 



-——-iujLq. (16) 



by tactically tuning the driving voltages. For instance, 
introducing asymmetry to the Gaussian pulse [i.e., set- 
ting K 7^ 1 in Eq. pS]) ] magnifies the values of I{uj) as 
cij ^ 0+, and thus reduces its relative errors. However, 
no matter what type of external voltage is applied, the 
HEOM approach should yield exactly the same G{uj) in 
linear response regime. This has been proved via exten- 
sive tests (see Sec. lIVI for the results with delta-function 
voltages) . 

From Fig. 0] it is revealed that at a finite temperature, 
Rq is no longer a universal constant, but depends on the 
system-bath coupling strength F. In particular, Rq devi- 
ates significantly from the "universal" value ro/2 when F 
is minute. To understand this temperature dependence of 
charge relaxation resistance, we plot in the inset of Fig. [4] 
the calculated Rq as a function of a dimensionless quan- 
tity T/T under various temperatures. All the curves are 
found to overlap each other, which indicates a general 
trend of coherent electronic dynamics. This thus reveals 
that the response current spectrum depends parametri- 
cally on the ratio T/T. 

As shown in the inset, Rq becomes drastically larger 
than ro/2 as F ^ 5T. This can be rationalized as follows. 
The QD level (eo = 0) possesses an intrinsic broadening 
with the magnitude F due to the dissipative interaction 
with the lead. At zero temperature, all the tunneling 
electrons are injected along the Fermi surface of the lead 
(/i = 0), and hence the open QD system is in full res- 
onance resulting in Rq = ro/2. At a finite tempera- 
ture, some electrons in the lead are thermally excited. 
However, as long as the energies of the excited electrons 
(or holes) remain in the linewidth of the QD level, i.e., 
(— F,F), the electronic dynamic coherence is still con- 
served. If the temperature is increased further so that a 
significant portion of thermally excited electrons cannot 
be covered by the aforementioned energy window asso- 
ciated with the QD level, the resonance is partially lost. 
Therefore, the mismatch between the energy of tunnel- 
ing electrons and that of QD level disfavors the electron 
transport process and gives rise to a much larger Rq. 
This is consistent with Fig. [5] where Re[Z(a;)] versus the 
dimensionless quantity lu/T for different values of F are 
plotted. It is shown that in the low frequency {uj < F) 
regime, Re[Z(w)] is rather insensitive to to (with the max- 
imal deviation of roughly 3% at a; = F) for all simulated 
cases. However, Re[Z(w)] assumes a larger average value 
as F becomes weaker. Therefore, we conclude that at a 
finite temperature, the charge relaxation resistance Rq is 
no longer a universal constant, but depends on both the 
temperature and the coupling strength F. 



Although Z{u!) diverges at uj — [due to the fact that 
I{io) = at w = 0], Re[Z(w)] = Rq is a. finite con- 
stant and independent of w, provided that w > 0+. The 
charge relaxation resistance Rq is deduced to be ro/2 [cf. 
Eq. ^], where tq — h/e^ is the resistance quantum. At 
finite temperatures the complex impedance of the open 
QD system is calculated by the HEOM approach together 
with Eq. (O. The resulting Re[Z(w)] at w ^ 0+ are plot- 
ted in Fig.[5]as a function of F. To ensure the accuracy of 
Z{uj), it is vital to have high precision for the calculated 
/(w), especially in the low to range. This can be achieved 



C. Off— resonant tunneling cases 

We now turn to the cases where eo 7^ ijl. It is important 
to point out that as eo rises away from jjl, a much larger 
number of exponential functions needed to be adopted 
for the expansion of C^„j,g (i) to guarantee the accuracy 
of the outcomes. It is inferred that the short-memory 
components of bath correlation functions play significant 
roles in the off-resonant cases. Figure [5] depicts the tran- 
sient current responses to an asymmetric Gaussian volt- 
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FIG. 6; Transient current responses to an asymmetric Gaus- 
sian voltage pulse (see the inset) for various off-resonant 
cases. The parameters are k = 4, Tc = 100 ps, Th — 22 ps, 
and the rests in unit of mcV: A = 10"*, T = 0.104, T = 0.1 
and W = 10. 



age pulse for various dot energies eo . The current profiles 
are analogous to the resonant cases shown in Fig. [21 ex- 
cept that the peak amplitude of /(t) declines drastically 
as Co deviates continually from /i. The corresponding 
frequency-dependent admittances are depicted in Fig. [71 
where both the real and inraginary parts of G{uj) ex- 
hibit conspicuous blue shifts with the increasing eg. Intu- 
itively, as the dot level goes up, fewer electrons can tunnel 
through the potential barrier at the system-bath inter- 
face and dwell on the QD, and hence the electrochemical 
capacitance of the reduced system becomes smaller. This 
coincides with the tendency shown in Fig.[7Ja), where the 
descending slope of Iin[G'(w)] in the low lo range becomes 
less steep as eo increases. 



IV. TRANSIENT ELECTRONIC DYNAMICS IN 
NONLINEAR RESPONSE REGIME 



As the applied voltage intensifies, the dynamic re- 
sponse of the open QD system goes beyond the linear re- 
sponse regime. In such a case, the associated frequency- 
dependent admittance depends explicitly on the specific 
type of the external voltage. It is thus difficult to derive a 
general expression for G{uj), and the electronic dynamics 
needs to be studied case by case. 

Since the HEOM approach admits arbitrary time- 
dependent applied voltages, no extra effort is required 
for calculations under higher external biases, i.e., the nu- 
merical procedures of Eq. ^ can be extended straightfor- 
wardly to the nonlinear response regime. In this section, 
transient currents driven by three types of turn-on volt- 
ages will be presented in the subsections: (a) an asym- 
metric Gaussian function, (b) a step function, and (c) a 
delta function. 
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FIG. 7: (a) Real and (b) imaginary parts of G{uj) for various 
off-resonant cases. Same parameters are adopted as in Fig. [6] 
In both (a) and (b), the lines represent different eo in unit of 
meV. 



A. Frequency— dependent admittance under 
asymmetric Gaussian voltage pulses 

In Fig.[51we plot the transient currents under asymmet- 
ric Gaussian voltages of the amplitudes A ranging from 
0.1 to 1 meV (recall that in linear response cases, A is set 
lower than 10~^ meV). It is observed that as A increases, 
the electron accumulation and depletion periods become 
separated from each other. Especially for A = 1 meV, the 
transient current almost vanishes in the time interval of 
150 ~ 200 ps, whereas in the due course A(i) still keeps 
rising. The corresponding frequency-dependent admit- 
tances are shown in Fig.[9l where two adjacent lines are 
separated vertically by 0.5 50 for clarity. It is evidently 
indicated that as the applied voltage increases, higher- 
energy current components are gradually activated, and 
the resulting G{uj) appears more fluctuating. The com- 
plicated lineshape of G(w) seems to exclude any simple 
equivalent classical circuit that can describe quantita- 
tively the coherent electronic dynamics of the open QD 
system. 



B. Transient current driven by a step function 
voltage and its spectrum analysis 

The exact time-dependent current I{t) driven by a 
step function voltage has been obtained by the NEGF 
method. ^^ Especially with the WBL approximation 
{W -^ cxd), its Fourier transform I{lu) can be evaluated 
conveniently via an EOM for the reduced single-electron 
density matrix for the reduced system32. Under a step 
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FIG. 8: Scaled transient currents under asymmetric Gaussian 
voltages of different amplitudes. The parameters are k — 4, 
Tc = 200 ps, Tj, = 44 ps, and the rests in unit of meV: T = 
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FIG. 10: Transient current responses to a step function volt- 
age. The lines represent different bandwidths. Other param- 
eters (in unit of meV): A = 5, T = 0.02, F = 0.05 and eo = 3. 
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function voltage, the lead level shift is A(i) — A8(i). 
With a tiny A, the linear-response admittance as well as 
the associated charge relaxation resistance Rq should be 
reproduced from the above derivations for I(uj)^ 

The transient currents calculated by the HEOM ap- 
proach are plotted in Fig.[Tni It is shown that with a 
narrower lead conduction band, the initial overshooting 
of I{t) becomes less prominent and it takes longer time 
for the open system to reach the steady state. High- 
frequency oscillations are clearly observed for I{t) at the 
time < t < 25 ps, which are due to the large amplitude 
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of applied voltage. This nonlinear effect is confirmed by 
making a comparison to Fig.[lja), where the rapid os- 
cillation is absent from the transient current due to the 
much smaller A. It is revealed from the Fourier analysis 
that the characteristic oscillation frequency is centered at 
ujQ = A — Co, a.s depicted in Fig.[TTJ In Fig. [T^ we further 
investigate I{lu) by plotting its dependence on eo under 
a fixed A. For all cases /(w) either reaches an extreme 
point or undergoes a sudden change in terms of its value 
at w = Wo- It is interesting to note that for eq < fj, the 
plotted Re[/(a;)] exhibits a peak at a; = loq, while for 
eq > /i a dip shows up. As eo is drawn closer to /i, the 
current response at the frequency loq is more accentuated. 
In the time domain, this implies an enhanced oscillation 
amplitude for the transient current. For the resonant 
case, i.e., eo = /i, the function form of Rc[/(a;)] near luq 
is ln(a;^ + ^)/^ with x — 2{uj — ujo)/T, which is different 
from the typical Fano line shape {x + q)"^ /(x'^ + 1),M 

The thermal influence on the transient current is also 
explored. In Fig. [13] we plot Re[/(a;o)] as a function of 
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temperature. It is found that the system under investi- 
gation starts to respond sensitively to the environmen- 
tal temperature at T ~ 0.01 meV. Upon further increase 
of T, the thermal effect overwhelms the system-bath- 
coupling-induced linewidth, which is 0.05 meV in this 
specific case, and thus dominates the electronic dynamics 
of the reduced system. 



C. Transient current driven by a delta function 
voltage and its spectrum analysis 

A delta function external field has been used to study 
the transient electronic dynamics of an isolated molec- 
ular system<2^ For an open system, such as the single- 
level QD system of our primary interest, a delta function 
applied voltage is also useful to investigate its dynamic 
properties. The time-dependent level shift considered 
is A(i) — AQ{t)S{t). The advantage of such a delta 
function voltage is two-fold: (1) its Fourier transform 
A(w) = A/2 is a constant, thus the current response 
of any frequency can be detected, with I{u!) oc G{llj); (2) 
the HEOM (fT2|) becomes a set of time-independent linear 
equations at t > 0+, and an efficient Chebyshev propaga- 
tor can be employed to solve the evolution of p„[t)r^^ 
which greatly reduces the computational cost. 

To simulate the transient current in response to a delta 
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function voltage by the HEOM approach, we notice that 
A(i) takes effect only within the infinitesimal interval, 
i.e., t G (0,0+). Therefore, the reduced dynamics at 
(t > 0+) can be solved by normal propagation of Eq. (fT^ 
(in absence of bias voltages) with the initial condition 
adjusted to p„{0~^). Consider a moderate amphtude A, 
so that all the ADOs remain finite within t G (0,0+). 
Therefore, neither /5n+* nor pj,~^ contributes as Eq. ([1 
is formally integrated from time i = to 0+: 

.0+ 

Pn(0+)=Pn(0)-fzVCTn / A(r)pn(T)dT 

"'0 



1-iA 



— Pn(0). 



(17) 



Here Pn(0) are the equilibrium reduced density matrix 
and associated ADOs, and A„ is expressed as follows. 



/•0+ ^ 

A, =^a,y A(T)dT--^(7n. 



(18) 



Due to the fact that An|fi=o = we have p{Q^) — p(0), 
i.e., the reduced density matrix remains continuous upon 
the delta function perturbation. 

In Fig.[T3]we plot the calculated transient currents cor- 
responding to different bandwidths. In both cases the 
current is instantaneously switched on to its maximal 
value at t = 0+, and then relaxes back to zero. The peak 
value Imax = -^(0+) is 21.4nA for W = 15meV, and 
11.5 nA ioi W — I meV, respectively. The corresponding 
current spectrums are depicted in Fig. 1151 The similarity 
in hneshape between Fig. [3] and Fig.[Tl]is noted, although 
the former is a linear response result while the latter be- 
longs to the nonlinear regime. This is actually due to the 
factorization property of I{uj) for the resonant tunneling 
case?^ where the current spectrum can be expressed as 
I{lu) = sin(A/2)X(ci;) in the resonant case (eg = /i = 0), 
where Xltu) is some complex function independent of A. 
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response to an asymmetric Gaussian voltage pulse, for 
both resonant and ofF-resonant tunneling cases. It is 
verified that at a finite temperature, the dynamic admit- 
tance of the open QD system in the low frequency range 
can still be characterized by a classical RLC circuit. 
However, the charge relaxation resistance Rq is found 
deviated from half a resistance quantum, but depends 
on both the temperature and the system-bath coupling 
strength. The concept of the equivalent classical circuit 
breaks down under higher bias. Complicated nonlinear 
features are observed, such as the activation of high- 
frequency current components with an increasing voltage 
amplitude. Transient current responses to step and delta 
function voltages are also explored. The basic features 
of the associated electronic dynamics are analyzed and 
discussed. The analytical and numerical results obtained 
in Sec, mil and llVI serve as convenient basis to understand 
the electronic dynamics of an interacting open system, in 
which the HEOM formalism is in principle exact. Work 
along this direction is underway. 

The HEOM approach also has great implica- 
tions for time-dependent quantum transport in 

nanoscopic molecular devices. Formalisms based 

on NEGF method ^1^8.29^.31^^^^^^^^ pj^, 

quet theoryfSa^^i^iMiis and QDTi&>ilii^>^ have been 

proposed. First-principles calculations have been carried 

out on realistic electronic devices 32^ However, in these 

simulations, approximate schemes for the dissipative 

dynamics in real time are inevitably introduced. For 

instance, the complete second-order formulation—!^^ 

and the WBL approximation were adopted to simulate 

the steady and transient current through a molecular 

device with time-dependent density-functional theoryi^ 

It is thus important that these approximate schemes 

can be improved systematically. Since the QDT-HEOM 

approach is capable of yielding exact results for quantum 

dissipative dynamics, it can be utihzed to cahbrate the 

approximated methodologies, and we also expect it to 

provide some guidelines for the potential progress of the 

approximate formalisms. 
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FIG. 15: (a) Real and (b) imaginary parts of /(oj) under a 
delta function voltage. Results with finite W correspond to 
7(t) in Fig. 1141 The common parameters are (in unit of meV): 
A = 1, r = 0.1 and eo = 0. 



V. CONCLUSIONS AND COMMENTS 

To conclude, we investigate the quantum coherent elec- 
tronic dynamics of a single-level noninteracting QD cou- 
pled to one electrode. Simulations are carried out based 
on the HEOM formalism of QDT.— In the linear response 
regime, quantitatively accurate frequency-dependent ad- 
mittance is obtained by the calculated transient current 



APPENDIX: NUMERICAL VALIDATION 

We verify our numerical implementation of HEOM for- 
malism by comparisons to known exact quantum trans- 
port results achievable via other methods. Three cases of 
noninteracting QDs coupled to left and right electrodes 
(L and R) of Lorentzian spectral density functions are 
demonstrated as follows. 

The first case studies the time-dependent transport 
through a single-level spinless QD driven by a step- 
function voltage pulse, calculated before exactly by the 
NEGF method,2ii^iS In Fig.[ll] we plot the transient 
currents calculated via the HEOM approach on the same 
system as that of Fig. 2 in Ref.i21,. Our results agree 
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FIG. 16: Transient current through R-lead, /^(f), in re- 
sponse to a step-function voltage pulse applied on R-lead. 
The lines correspond to different lead bandwidths. Other pa- 
rameters: Tl = Tfl = 0.5 r, r = 0.1 r, Wl^Wr^W and 
A = lOr. This figure reproduces Fig. 2 in Ref.[2l|. 
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FIG. 17: Transient currents through L-lead, I hit), in re- 
sponse to a step-function voltage pulse applied on R-lead. 
The lines correspond to different voltage amplitudes. Other 
parameters: Tl = Fa = 0.5 F, T = 0.05 F and Wl ^ Wr ^ 
20 F. This figure reproduces Fig. 1 in Ref. ISa. 
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FIG. 18: Time-averaged current / (in unit of /max = eV /2h) 
through a noninteracting double QD, driven by ac gate volt- 
age (see Ref. l40| for details) . The interdot coupling strength 
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resonance energies lon ~ -^ 1/(^2 — ei)^ + 4k;2 are labeled b 
numbers A'^ = 1 ^ 4. This figure reproduces Fig. 3 in Ref. [4] 
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FIG. 19: Transient current /^(i) (in unit of /max ~ eF/2ft) 
through a noninteracting double QD. The system setup and 
other parameters are same as in Fig. 1181 



quantitatively. The QD is initially in equilibrium with 
zero bias. External voltage pulses are switched on from 
the time t = 0, which results in current flows through the 
leads L and R a,t t > 0. Hereafter we denote ^L{t) and 
Aii{t) as the time-dependent energy shifts for the lead 
L and R due to the applied voltages Vl(*) and ^^(i), 
respectively, i.e., Aa(t) — —eVa(t). For simulations pre- 
sented in this section, we set ^Lit) — for all t, and 
^ii(*) = ^0(0: where Q{t) is a step function turned 
on at t = 0, and A is a constant amplitude. The system 
level energy is time-dependent as eo{t) — eo + A£)(t) with 
Az3(t) = [AL(t) + Afl,(i)]/2. 

The second example is similar to the first one, but 
with A£i(t) = 0. The calculated transient currents dis- 
played in Fig.[T7] accurately reproduce previous simula- 
tions (c/. Fig. 1 in Ref.'52') by the NEGF method with all 
the non-Markovian features preservedi^^iS 



The third example investigates the resonant photon- 
assisted tunneling through QDs where results have been 
obtained by a combined method of NEGF and Flo- 
quet formalisms. ^"'^^ In Ref. 40. nonequilibrium electron 
pumping through a double QD system driven by ac gate 
voltage was simulated. Multi-photon-assisted tunneling 
was resolved in time-average current spectrum. Since the 
HEOM approach admits an arbitrary time-dependent 
external field, a sinusoidal gate voltage can be treated 
readily. However, instead of the frequency-domain cal- 
culation conducted in Ref.kol. with our present coding 
scheme we need to propagate Eq. p^ in real time for 
every individual frequency. This makes the simulation 
for the entire current spectrum rather tedious. Nonethe- 
less, to further justify our numerical procedures, we have 
managed to carry out such a simulation and compared 
to reported results. A noninteracting double QD sys- 
tem coupled to two leads is studied. The time-averaged 
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current is evaluated as the system reaches a quasi-steady 
state for each gate voltage frequency uj. We take the same 
parameter set as adopted by Fig. 3 in Ref.[40, and the 
HEOM calculation result is plotted in Fig.[THl The quan- 
titative agreement between our Fig. [18] and result shown 
in Ref.|40| is noted. In particular, the characteristic A^- 
photon-assisted resonance frequencies are correctly re- 
produced, as demonstrated in Fig.[THl This test again 
validates our numerical procedures. It is worth mention- 



ing that the Floquet formalism treats the quasi-steady 
dynamics driven by an ac external field, while our HEOM 
formalism allows access to much broader information be- 
yond this, such as the establishment of a quasi-steady 
state in real time. For instance, Fig. flW a) and (b) depict 
the time evolution of a double QD driven by an ac gate 
voltage switched on from t = for uj/w = 0.1 and 1.1, 
respectively, where w is the inter-dot coupling strength 
as the system setup is the same as in Fig. [TBI 
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